General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



I 


(IASA-CB-176035) IVTEiACTJVA A1BCBAPT DCS- ‘1064 

PLIGHT COITIOL AMD ABBOELAS1IC SI M)ILI2ATICN 
Semiannual fleport, 1 Mov. 1964 - 30 Apr. 

1905 (Purdue Oniv.) 20 p bc A02/HP A01 Onclas 

CSC I 01C Ci3/0 8 21805 

INTERACTIVE AIRCRAFT FLIGHT CONTROL AND 

AEROELASTIC STABILIZATION 


NASA/Langley Research Center Grant-NAG-1-157 
Semi-annual Report 

1 November 1984 through 30 April 1985 


Submitted by: 

Dr. Terrence A. Weisshaar 
Principal Investigator 


SCHOOL OF AERONAUTICS AND ASTRONAUTICS 
PURDUE UNIVERSITY 
WEST LAFAYETTE, INDIANA 47907 



Summary 

This report covers research performed under sponsorship of 
NASA/Langley Grant NAG-1-157 during the period 1 November 1984 through 
30 April 1985. An analytical model of a 3-D airfoil was used to study 
an optimization procedure formulated to enhance stability of an airfoil 
through Integrated structural and control synthesis. This procedure is 
discussed In this report, together with preliminary results. These 
results show that a sensitivity derivative approach utilizing structural 
parameters, weighting matrix parameters and optimal control parameters 
(In this case, the design airspeed) is effective in determining the 
"best" structural/control design. 

Two trips were taken during this reporting period. The first trip 
was taken 28-29 January 1985 to visit Langley Research Center to present 
a research progress report. The second trip was taken from 1 April to 3 
April 1985 to attend the Second International Conference of Aeroelastl- 
city and Structural Dynamics In Aachen, West Germany. A paper was 
presented at this conference. A trip to MBB Aircraft In Munich was also 
Included from 4 April to 5 April 1985 to discuss common research 
efforts. A trip summary is Included In this report. 


Discussion 

The objective of the current work was to develop and to demonstrate 
a procedure to ensure aeroelastlc stability of an airfoil up to and 
Including a target airspeed. To accomplish a portion of this objective, 
optimal control techniques were used to design the active control sys- 
tem. In addition, sensitivity derivatives were computed to assess how 
changes In system parameters affect the optimal control law design. In 
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addition to control parameters auch as the elements of the state weight- 
ing matrix (the [Q] matrix) and the design speed at which the control 
lav is formulated, denoted aa , the position of the shear center, 
expressed as a parameter a , was used as a design variable. U is a non- 
dimensional airspeed, while a^ is a ncnd linens lonal coordinate described 
below. 

The objective of the new procedure Is to define a control law such 
that the aeroelastlc system is stable at all speeds below a certain 
airspeed, called U g . A problem that arises with the use of optimal con- 
trol techniques to define such a control law Is that the system may be 
stabli! at the airspeed at which the control law Is formulated, but may 
be unstable at lower airspeeds. As a result the design airspeed for the 
control law does not necessarily correspond to the maximum airspeed to 
be reached before aeroelastlc stability is encountered. Conversely, 

U may be used as a parameter in the stability augmentation problem. 

Similarly, the elements of the state weighting matrix [Q^l used 
for optimal control synthesis are arbitrary. However, the choice of 
these elements affects the control law and the off-design perfor- 
mance of the active control. Setting these parameters in an "optimal" 
manner is also advantageous. Therefore, the elements become duslgn 
parameters. 

The structural design affects the behavior of the flexible struc- 
ture. The example chosen for the current research is a 3-degree of 
freedom, typical section model as shown in Figure 1. The degrees of 
freedom are airfoil pitch, a, airfoil plunge, h, (bending) and control 
surface rotation, 0. In particular, for this 3-DOF model, the position 
of the shear center affects the aeroelastlc stability of the airfoil. 
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Figure 1 - Three-degree- of-freedom typical section airfoil. 

In reality, the airfoil skin thicknesses, span areas and locations all 
affect the shear center location. A parameter, a g , Is used to denote 
shear center position. The parameter a^ is the nondlmenslonal position 
(position divided by semi-chord) of the airfoil shear center with 
respect to the airfoil midchord. As a result, a g ■ 0 represents the 
shear center position at the airfoil midchord, while a g - -0.50 
corresponds to a shear center placement at the airfoil quarter chord. 

For this study we have assumed that there is no weight variation 
(penalty) to be Incurred when the shear center Is moved. Furthermore, 
the sectional center of mass Is held fixed as a^ is changed. These res- 
trictions do not Invalidate the results to be discussed nor th» pro- 
cedure used to generate t tese results. 

The shear center location has a substantial effect upon the "open- 
loop" aeroelastlc stability of the 3-DOF airfoil. The shear center 
position also has a substantial effect upon the value of the "cost func- 
tion" used to generate control laws using full-state feedback with 
optimal, steady-state, linear-quadratic-regulator theory (LQR theory). 
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Figure 2 - Constant optimal control cost contours for an actively controlled 

3-DOF typical section as a function of U and a . 

Des e 

In Figure 2 the effect of the parameter a g upon the standard LQR 

coat function J for different values of U_ is indicated. The open- 

Des 

loop flutter and divergence speeds, as functions of a g , are indicated as 
bold solid and dashed lines, respectively, in Figure 2. 

If a^ - -0.20 and Up eg - 6.0 an optimal control law may be designed 
such that closed-loop stability is ensured at U - 6.0. As seen in Fig- 
ure 2, the cost, J, is approximately J ■ 10. If Up es Is fixed at 6.0, 
but a e is changed to equal -0.4 (a forward shear center movement), sta- 
bility of the closed-loop system is still ensured at U ■ 6.0. However, 
the "cost" has increased to J ■ 20. While the control cost can not be 
directly translated into real costs, the implication here is that the 
system requires more effort to control. This is due to the close prox- 
imity of the open-loop stability boundary when a^ ■ -0.40. 

Next, let us examine the closed-loop behavior of this system. To 
generate Figure 3, active control laws were formulated, using LQR theory 
with U fixed at 6.0, and at several different values of a . Both 
cpen-loop and closed-loop instability speeds are shown in Figure 3. 
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Figure 3 - Closed-loop stability boundaries for an actively controlled 

3-DOF typical section. The design airspeed Is - 6.0 while 

control laws are formulated as various values of a 8 . 

e 

With fixed at 6.0 the closed-loop stability behavior may be 

degraded for one of two reasons. First of all, the open-loop stability 

- boundary may be degraded at speeds above Ujj e8 ■ 6.0 when active control 

is added. This is seen to be the case when a g is in the vicinity of 

-0.30. In a second case, when a is less than -0.40, a low-speed insta- 

e 

bility region appears due to the addition of feedback control. 

Notice also, that if one were to fix and maximize the airspeed 

at which closed-loop aeroelastlc instability occurred, this maximum 
would occur near a fi ■ -0.20. This value is slightly less than the max- 
imum obtained using passive means (change in a g only). 

In summary, the choice of both a g and Up es has a substantial effect 
upon closed-loop system stability. In addition, low values of the LQR 
coat function, J, do not necessarily lead to acceptable dynamic response 
at all speeds below the instability speed. As a result of these obser- 
vations a modified procedure was investigated to both take advantage of 
LQR theory, optimization theory and optimal sensitivity derivatives. 

This method was proposed by Mr. T.A. Zeiler who has included work by Mr. 
M.G. Gilbert to produce a hybrid scheme for selecting the best 
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structural/control design 


An Integrated Design Synthesis Scheme 
Let us define as our design objective the Increase in the aeroelas- 
tlc instability speed (either flutter or divergence) to a certain value, 
U g . This Increase will be accomplished by utilizing two sets of design 
variables, those involving structural parameters (in the present case, 
only a g ) and those involving active control law parameters. The optimal 
control design airspeed Up eg does not necessarily correspond to the 
design objective airspeed, U g . We also require closed-loop system sta- 
bility at all airspeeds below U g . 

At U or at any other airspeed below U , the closed-loop system 
8 8 

will yield eigenvalues of the form 

X t - + jw i (1) 

Our design objective requires that all values of o^ be negative (in the 

left half-plane) at speeds below U . A function, F , may be defined 

8 8 

such that 

1 N p0 i 

F -i ln[ I e *] (2) 

8 p i-1 

F g is called a cumulative constraint function. In Eqn. 2, p is a con- 
stant chosen to scale the problem while N corresponds to the number of 
distinct values of obtained from the closed-loop eigenvalue analysis. 
It can be shown (see Sobleskl, et al. , AIAA Paper No. 83-C832-CP) that 
F g is bounded as follows: 

o^(max) < F g < o^(max) + InN (3) 

Thus if at a given airspeed, U, F g is negative, then o^(max), the larg- 
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est value of at thla airspeed, U, la also negative and the system la 
a table. 

What happens If, at or some other airspeed U^, the value of F # 
la found to be positive? How do we modify the system control law and 
structure to reduce the value of F g to an acceptable value In an 
"optimal" manner? Certainly we could rely exclusively upon active con- 
trol. However, situations such as shown in Figure 3 might arise, for 

which "sub-crltlcal" instabilities appear. 

% 

Let us designate F g as being a function of a set of "design" param- 
eters, p. , in this case a , U_ and 0 . The functional F Is also a 

* r i* e Des 1J s 

function of the control, u. 

First, we select a finite set of airspeeds at Intervals up to U^; 
call these airspeeds U^. These speeds are chosen to monitor system sta- 
bility at "sub-crltlcal" airspeeds. At each airspeed except U , the 

8 

actively controlled closed-loop eigenvalues X^ yield functionals F^ such 
that 


Fj < 0 at Uj (4) 

where the form of F^ is the same as that Indicated in Eqn. 2. At our 
aeroelastlc design airspeed U^, F g > 0. We wish to choose a new set of 
design parameters, p^. In such a way that we minimize F g while keeping 
Fj < 0. In addition we require that the LQR cost function J be a 
minimum with respect to a choice of control, u at Changes In the 

parameters p^ to reduce F # so as to stabilize the system at will, In 
general, change the values of F^ at the intermediate design airspeeds. 
This will occur because of reshaping of the root locus curves for the 
closed-loop system. 
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The procedure Co be used to optimize the actively controlled system 
Is a modification of the multi-level, linear .decomposition procedure 
suggested by Sobleskl and co-workers. In the present case the system Is 
decomposed Into a structural subsystem and an active control subsystem. 
Let us first discuss the active control subsystem. 

At the active control level, the subsystem optimisation Is the 
solution to the optimal steady-state LQR problem, 


with 


min J ■ / [X C QCX + u Ru]dt 
u o 


(5) 


X - AX + Bu , 


( 6 ) 


The solution Is: 


u - GX - — R~ 1 B*PX , (7) 

where P Is the solution to the steady-state matrix Rlccatl equation. 

The change In J with respect to design parameters Is achieved through 
differentiation of the necessary conditions satisfied at the minimum J, 
known generally as Kuhn-Tucker conditions. The results of Interest are 
derivatives of the optimal feedback gain matrix wl'h respect to the 
design parameters, (assuming B * B(p)) 


3G - an* 3P /at 

9 ? R ® ^ • (8) 

That Is, the derivatives are constrained to describe how the optimal 
gains change with changes In parameters. Thus, the control-augmented 
state matrix (denoted as A + ) derivatives are written as: 


9A + 3A 3G 

dp dp dp • 

Eigenvalue sensitivity derivatives are found from 


( 9 ) 


\ 


8 


/ 

ax _ 3A 

1^" DIAC[E W E] , (10) 

where, X^ ■ diagonal matrix of system eigenvalues; E - eigenvector 
matrix of A + . 


These derivatives describe how the eigenvalues of the optimally con- 
trolled system change with changes In the design parameters. The real 
parts of the eigenvalues and their sensitivity derivatives, 


ao ax 

IT " <1T> • 


(ii) 


are. then used In the system level optimization. At this point we can 
change parameters such that we still have an optimally controlled 
system, that Is, at a minimum control "cost". 

For Lhis example there Is no structural cost, that is, no addi- 
tional weight is associated with changes In a^. 

At the system level, the optimization problem can be written in 
linearized form as: 


subject to 


min F 
*P« ' 


F + 
80 


m 

£ 

1-1 



( 12 ) 


m 3F 

(a) F - F + £ Ap < 0 (13) 

J J o 1-1 3p l 1 

(b) J - min J 

u 

(c) x - Ax + Bu 

where m - the number of active design parameters. Note that constraints 
(b) and (c), representing the subsystem optimization, are implicitly 
satisfied since optimal sensitivity derivatives of the eigenvalue real 
parts are used In construction of the Taylor series approximations to 
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the cost function F # and constraints. The derivatives In Eqns. 12 and 

13a ara given aa: 


N po. Jo. 

». 1-1 * ^ 


3p 


N po. 
E e 
1-1 


(U) 


To test this procedure a ’’simplex" algorithm adapted from linear 
optimization was used. The parameters - a^ and p^ - Up e8 were chosen 
as design parameters. An Initial design with a - -0.4 and U - 6.0 
was chosen. In this case the closed-loop system was unstable at U - 
7.0. The objecclve was to stabilize the system at - 8.0 In an 
optimal manner. To do this, the functional associated with ’J ^ - 7.0 was 
first reduced. Then the functional associated with - 8.0 became 
and It was reduced. 

Figure 4 shows the history of a^ and Up e8 versus number of design 
cycles, together with the values of F # and F^. After seven Iterations, 
no meaningful stability Improvements were possible. Figure 5 shows the 
system root locus for the final design. The final design appears to be 
a compromise between flutter of the closed-loop system in two different 
modes, that Is, a cusp In the flutter boundaries such as shown In Figure 


Figure 6 shows a plot of pole-zero "sdgratlons" with airspeed for a 

fixed value a - -0.4. At this value of a there Is close pole-zero 
e e r 

proximity (not exact cancellation) for an unstable mode. This may 
account for the fact that the control designed at this point does not 
produce superior stability boundaries. However, the process of reducing 
F ( moves the system away from this divergence region. 
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Figure 



Numier of Design Cycles 


4 - Values of s end U versus design cycle iteration 

number. AlSo shown Ire the values of the cumulative constraint 
functions F and F for ±ach design cycle. Integer numbers refer 
to associated values of U.. 
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Figure 5 - Root locus plot for final design cycle, a - -O.J_89 and 
U - 6.01S. Note that instability occurs at U ■ 8. AO 
an8*U - 8.43. 
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Figure 7 shows the pole-zero migration as a function of U for a 
configuration close to the final design. This figure shows that the 
problem of unstablllzablllty (uncontrollability of an unstable mode) is 
avoided by stabilization of the unstable mode through structural modifi- 
cation. There does appear to be some additional pole-zero separation, 
representing enhancement of controllability (as opposed to avoidance of 
Incontrollability) . 

In addition to the above study, an investigation was made of the 

effects of changing the diagonal elements of the output weighting 

matrix, Qjj* All previous studies had taken [Q^] to t> e an identity 

matrix. Figure 8 shows stability boundaries versus a g for * 6.0 

and airfoil pitch weighting of 0 - 100. With this weighting, the 

o 

closed-loop divergence region for large negative values of a g that 
exists when [Q^jl is an identity matrix is eliminated. (Compare Figure 
8 with Figure 3.) However, the flutter boundaries and the divergence 
boundary when a is near zero are unaffected by this change. The 
closed-loop divergence boundary appears to be associated with the con- 
trol deflection and, in fact, merges with the flutter boundary Just 
above it. Figure 9 shows stability boundaries for the control deflec- 
tion weighting, Q - 100, with all other diagonal elements equal to 
8 

unity. The divergence boundary for large negative a g Is again elim- 
inated while the main flutter boundary is only slightly affected. The 
divergence and flutter boundaries for a g near zero are shown to be more 
heavily Influenced by changes in Q . 

P 

Summary of Results and Future Work 
During the past six months a procedure has been developed to 
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Figure 7 - Plots of actively controlled system poles and zeroes as a function 
of airspeed for a 3-DOF model with a - -0.2. Real axis poles 
and zeroes are not shown. e 
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Figure 9 - Open-loop 
Q 6 - 100, 
system at 


and closed-loop stability boundaries with Q 

Q. . ■ 0, 1 * J for an actively controlled 

various values of a . 

e 
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redesign the structure end the control system to augment the aeroelsstlc 
stability of a idealized aeroelastic system- A number of problems 
remain to be resolved and some procedures need to be formalized. How- 
ever, it appears that a major step toward the integration of the 
structural/control optimization has been accomplished. 

Future efforts will be directed towards exercising the method 
further on the 3-DOF model as well as the 4-DOF model with "fuselage" 
pitch freedom. 

Work will also begin on more realistic models that incorporate 
multi-mode, laminated composite structures. A Master's degree student, 
Mr. V.J. Sallee will begin a 10 week residency at Langley Research 
Center in mid-May 1985 to learn to operate the ISAC code so that these 
new thrusts can begin. 


Trip Report 


Second International Conference on Aeroelastlclty and Structural Dynam- 
ics, Aachen, West Germany, 1 April through 3 April 1985. 

A trip was made to this conference to present a paper entitled 
"Tailoring for Aeroelastlc Stability and Lateral Control Enhancement". 
The conference was well-attended by European specialists In this area of 
research and technology. In the field of aeroelastlc tailoring the 
presentations were for the most part not state-of-the-art as we know It 
In the USA. Personal discussions with engineers from the German air- 
craft establishment lead one to believe that their capabilities and 
Interest are much greater than conference presentations would Indicate. 
Some mention was made of the aeroservoelastlc tailoring problem although 
Ideas about how to approach the subject were not forthcoming. The 
Israeli aircraft engineers and researchers seem also to have a submerged 
Interest In tailoring. 

A visit to MBB, Munich was made on April 4 and 5. Mr. Otto 
Sensburg was the host for this visit. During this visit there was a 
keen interest expressed by the Germans In Innovative technology. I was 
briefed on their design for a tailored composite vertical stabilizer and 
planned scale-model tests. They also plan active control testing. A 
lasting Impression was that of a highly trained and qualified group of 
German engineers with many new Ideas, not all of which they were willing 
to discuss. 
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